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OBJECTIVE 


We  evaluate  and  characterize  seismic  data  from  digital  network  and  array  sites  in  and 
around  the  Former  Soviet  Union,  in  order  to  provide  new  constraints  on  seismic  wave 
propagation  characteristics  in  the  Caucasus  Ranges,  Middle  East  and  Iranian  Plateau. 

Strong  changes  in  the  attenuation  and  velocity  behavior  of  seismic  waves  are  associated 
with  the  mountain  belts  here,  and  make  it  difficult  to  infer  wave  behavior  from  experiences 
gained  within  stable  cratons.  Seismograms  from  the  JSP-Caucasus  Seismic  Network 
collected  since  1991  include  a  large  number  of  local  and  regional  events  in  the  region. 

These  waveforms  are  used  for  assessing  the  ampUtude  and  phase  of  major  wave  groups 
along  regional  paths,  and  constraining  earth  structure  beneath  the  network.  Results  provide 
information  on  the  behavior  of  common  event  discriminants  and  detection  thresholds  for 
regional  paths. 

RESEARCH  ACCOMPLISHED 

Severe  changes  in  wave  behavior  are  observed  between  paths  that  cross  the  Caucasus 
and  other  Alpine-Himalayan  mountain  belts,  and  those  that  do  not  (Fi^re  1).  These 
changes  represent  substantial  attenuation  of  phases  often  used  to  discriminate  explosions 
from  earthquakes  and  to  estimate  event  size  and  type,  such  as  Lg  and  other  high-frequency 
surface-related  phases.  Work  to  date  has  concentrated  on  quantifying  the  behavior  of 
seismic  waves  as  they  cross  the  Caucasus,  mostly  by  examining  phase  amplitudes  for 
events  within  10°  of  the  Caucasus  Seismic  Network  (CNet).  Other  facilities,  such  as  the 
Geyokcha,  Turkmenistan  array  also  provide  digital  seismograms  that  can  give 
complementary  information. 

/.  Body  wave  attenuation  across  the  Greater  Caucasus  range 
front 


Summary.  A  regionally  abrupt  boundary  in  seismic  attenuation  is  documented 
between  the  Russian  platform  and  the  tectonically  active  Greater  Caucasus.  Attenuation 
( Q.  ' )  is  measured  from  body  waves  recorded  by  a  digital  seismic  network  that  straddled 
north  flank  of  the  Greater  Caucasus  from  1992  to  1994.  Spectral  fall-off  of  body  waves  is 
measured  from  186  events  within  4°  distance  of  the  stations,  and  a  frequency-independent 
attenuation  parameter  ( /* )  and  source  comer  frequency  is  determined  from  fits  to  amplitude 
spectra.  Most  fits  are  at  frequencies  between  2  and  15  Hz,  limited  by  background  noise. 
These  t*  measurements  are  then  inverted  for  Q~' ,  for  parameterizations  of  varying 
complexity.  Results  show  (2j  of  775  ±75  and  2060  ±315  for  mountain  and  platform  paths, 
respectively,  and  Qp  ~  Qs-  Inversions  show  more  than  30%  variance  reduction  for 
parameterizations  that  allow  (2~‘  to  vary  laterally  relative  to  laterally  homogeneous 
parameterizations.  Site  corrections  also  significantly  reduce  variance  of  fit  between 
predicted  and  measured  t*.  High  attenuation  is  observed  beneath  high  topography, 
separated  by  a  sharp  boundary  from  the  low-Q"'  platform.  Comparison  of  body  wave 
spectra  to  those  of  coda  suggests  intrinsic  absorption  rather  than  scattering  controls  the  Q 
measurements,  at  least  beneath  the  mountains.  If  the  observed  Q~'  contrasts  are  due  to 
changes  in  temperature,  they  suggest  70-400°  C  higher  temperatures  beneath  the  mountain 
relative  to  the  platform.  These  temperature  increases  may  not  be  high  enough  to  generate 
widespread  partial  melting  beneath  the  mountains,  but  are  sufficient  to  produce  regional 
metamorphism,  and  can  contribute  substantially  to  isostatic  compensation  of  the  mountains. 
Whatever  the  cause,  these  data  show  that  the  boundary  in  seismic  attenuation  is  abrupt  and 


large  between  stable  shield  and  an  adjacent  mountain  belt,  showing  that  Q  is  a  sensitive 
indicator  of  tectonic  process. 

I.l  Introduction 

Eurasia  is  dominated  by  a  first-order  distinction  between  stable  shields  and  active 
mountain  belts  associated  with  Alpine-Himalayan  collision.  The  tectonic  style  of  these 
mountain  belts  vary  considerably  (e.g.  Jackson  &  McKenzie  1984),  but  in  many  cases  the 
boundary  between  stable  Eurasian  cratonic  regions  and  the  mountains  to  their  south  are 
abrupt.  For  example,  peaks  of  the  Greater  Caucasus  exceed  5  km  in  elevation,  within  view 
of  the  Russian  (Scythian)  platform  where  elevations  are  uniformly  below  500  m  (Figure  1). 
At  long  wavelengths  (>1000  km),  it  has  long  been  known  that  differences  between  the 
deeper  structures  of  shields  and  mountains  are  also  large.  Seismic  wave  propagation  over 
regional  distances  beneath  stable  Eurasia  is  very  efficient  (Zhang  &  Lay  1994;  Mitchell,  Xie 
&  Pan  1996;  Kim  1987),  while  many  seismic  waves  transmit  poorly  through  adjacent 
mountains  (Kadinsky-Cade  et  al.  1981;  Rodgers,  Hearn  &  Ni  1997). 

It  is  less  clear  how  these  first-order  transitions  in  wave  propagation  manifest 
themselves  over  scales  of  tens  of  km,  scales  at  which  many  tectonic  processes  operate. 
Transitions  at  depth  may  correspond  to  geologic  boundaries  on  the  surface,  to  topographic 
boundaries,  or  neither.  For  example,  across  the  Rocky  Mountain  front  of  North  America, 
it  has  been  found  that  subsurface  boundaries  are  offset  by  hundreds  of  kilometers  from  the 
topographic  range  front  (e.g.  Lee  &  Grand  1996;  Humphreys  &  Dueker  1994).  Here,  we 
take  advantage  of  new  seismic  data  across  one  such  boundary,  that  between  the  Greater 
Caucasus  and  Russian  platform  to  explore  the  subsurface.  Because  the  effect  of  attenuation 
is  strong  (Figure  2)  and  well  sampled,  we  use  attenuation  variations  to  describe  the  shield- 
mountain  transition.  A  second  goal  of  this  study  is  to  better  understand  lateral  variations  of 
signal  propagation  characteristics,  which  play  a  key  role  in  seismic  monitoring  at  regional 
distances.  We  find  that  large-scale  boundaries  between  tectonic  provinces  can  generate 
sharp  boundaries  to  seismic  attenuation,  and  that  the  boundary  in  the  Caucasus  at  least 
correlates  well  with  that  seen  in  topography. 

The  method  used  here  to  determine  Q  is  based  on  measurements  of  the  spectral  fall-off 
of  body  waves,  similar  to  Hough  et  al.  (1988a).  This  approach  gets  around  many  of  the 
problems  associated  with  Q  determinations  based  on  absolute  amplitudes,  which  can  be 
greatly  affected  by  near-surface  effects,  focusing  and  defocusing,  and  source  behavior 
(e.g.  Frankel  &  Wennerberg  1987).  Body  waves  also  sample  energy  that  travels  close  to 
the  ray-theoretical  phase  arrival  time,  so  that  the  sampling  volume  is  small,  unlike  the  large 
volumes  sampled  by  coda  waves.  Hence  resolution  of  structures  is  much  greater  than  for 
coda  techniques,  and  we  can  resolve  a  sharp  attenuation  boundary  within  the  crust  that 
coincides  with  the  boundary  between  tectonically  active  and  inactive  regions.  The 
sharpness  of  this  boundary  suggests  that  long-wavelength  Q  measures  may  be  severely 
aliases  across  mountain  fronts.  As  well,  the  mechanisms  for  producing  low  Q  may  be 
closely  tied  to  the  mechanisms  for  producing  topography  in  collisional  belts. 

1.2  Regional  Setting  and  Seismicity 

The  Greater  Caucasus  accommodate  Cenozoic  convergence  between  Eurasia  and 
Arabia,  currently  29  mm/yr  at  N26“W  in  this  region  (DeMets  et  al.  1990).  The  mountains 
strike  SE-NW  between  the  Black  and  Caspian  seas,  and  form  a  part  of  the  continuous 
orogenic  belt  along  the  southern  edge  of  the  Russian  Platform.  The  southern  edge  of  the 


Platform  is  covered  with  Mesozoic  and  Cenozoic  sediments,  overlying  Hercynian  basement 
of  the  Scythian  craton  (Adamiya  et  al.  1992).  Uplift  and  convergence  began  in  the  Middle- 
to-Late  Miocene,  and  intensified  during  Mid-Pliocene  (Burtman  1989).  The  core  of  the 
Greater  Caucasus  consists  of  Precambrian  and  Paleozoic  igneous  and  metamorphic  rocks 
(e.g.  Zonenshain,  Kuzmin  &  Natapov  1990),  while  the  southern  slope  of  the  mountains  is 
a  thrust  belt  involving  backarc-type  sediments,  deposited  throughout  the  Jurassic, 
Cretaceous  and  Paleogene  (Zonenshain  &  Le  Pichon  1986).  North-dipping  thrust  sheets  or 
nappes  characterize  extensive  deformation  on  the  southern  slope  of  the  Greater  Caucasus, 
while  a  gently  north-dipping  monocline  bounds  the  northern  slope  (e.g.  Philip  et  al.  1989). 
T.atp  Quaternary  volcanism  occurs  in  isolated  locations  throughout  western  part  of  the 
Greater  Caucasus  (Philip  et  al.  1989),  including  Mt.  Elbruz  at  the  SE  edge  of  the  network 
(Figure  1),  which  may  be  still  active  (Kovalev  &  Masurenko  1972). 

The  region  is  seismically  active,  with  the  majority  of  seismicity  concentrated  to  the 
south  of  the  Greater  Caucasus  near  the  Iranian  plateau  (e.g.  Jackson  &  McKenzie  1984). 
The  April  29, 1991  Racha  earthquake  (moment  magnitude  Mw=7 3)  is  the  largest 
instrumentally  recorded  earthqu^e  in  the  Greater  Caucasus  (Triep  et  al.  1995;  Gorshkov 
1984).  The  focal  mechanisms  of  that  and  smaller  earthquakes  show  thrust  faults  striking 
parallel  to  the  mountains,  with  shallow-dipping  planes  dipping  to  the  north  (Jackson  & 
McKenzie  1984).  Thus,  the  north-dipping  thrust  faults  that  bound  the  southern  slopes  of 
the  Greater  Caucasus  are  still  active  (Triep  et  al.  1995).  Many  of  the  events  studied  here 
are  near  the  1991  aftershock  zone  (near  43“E,  42*N).  East  of  45°E  most  earthquakes 
studied  here  occur  north  of  the  mountains,  where  a  south-dipping  thrust  belt  is  exposed 
(Figure  1). 

The  Caucasus  are  a  region  of  anomalously  low  seismic  velocities  and  high  attenuation 
(low  Q).  Rautian,  Khalturin  &  Shengeliya  (1979)  show  coda  Q  value©=near  330  at  1  Hz. 
Kadinsky-Cade  et  al.  (1981)  and  Rodgers  et  al.  (1997)  both  show  that  at  regional  distances 
(5-25°)  the  phase  5„  propagates  weakly  or  not  at  all  at  frequencies  near  1  Hz,  through  the 
Caucasus.  As  well,  Lg  is  shown  in  these  studies  to  propagate  poorly  across  the  Black  Sea, 
the  southern  Caspian  Sea,  and  through  some  regions  of  high  topography.  The  strong 
attenuation  led  these  authors  to  postulate  the  existence  of  widespread  partial  melt  beneath 
the  Caucasus  within  the  uppermost  mantle.  A  similar  inference  was  made  by  Heam  &  Ni 
( 1994)  who  observed  extremely  low  P„  velocities,  near  7.6  km/s,  associated  with  the 
Caucasus.  Our  results  confirm  and  quantify  the  low  attenuation  values  associated  with  the 
Caucasus,  and  document  the  northern  boundary  of  the  attenuative  region.  We  show, 
however,  that  the  data  are  consistent  with  elevated  temperatures  and  do  not  require  that 
partial  melt  be  present. 

1.3  Data  and  Methodology 

Data  were  recorded  by  the  Greater  Caucasus  seismic  network,  operated  from  199 1 
to  1994  through  a  cooperative  project  among  IRIS,  Lamont-Doherty  Earth  Observatory  and 
the  OME-Obninsk  group  (Abers  1994).  Stations  were  located  on  the  northern  slope  of  the 
Greater  Caucasus,  and  straddled  the  transition  from  the  Greater  Caucasus  to  the  Russian 
Platform  (Figure  1).  The  northern  stations  (NAG,  MIC,  KUB)  lie  on  a  flat-lying,  poorly 
consolidated  Neogene  sediments,  while  southern  stations  (KIV,  GUM,  KNG)  lie  on  a 
tilted,  well-lithified  Mesozoic  sandstones  and  limestones.  The  network  consisted  of  up  to 
six  digitally  telemetered  three-component  intermediate-period  sensors  (Kinemetrics  SV/SH- 
1),  digitized  at  60  samples/s  (sps)  by  Nanometrics  RD-3  digitizers.  The  configuration  has  a 
flat  response  to  velocity  at  frequencies  between  0.2  and  24  Hz.  The  five  remote  stations 


digitally  transmitted  the  signals  via  a  one-way  VHF  telemetry  links  to  a  central  station 
Kislovodsk  (KIV),  where  network-triggered  data  streams  were  recorded  and  initially 
examined. 

A  total  of  1571  events  were  located  from  January  1992  to  first  half  of  1994,  of  which 
186  events  were  selected  for  study.  Most  of  these  186  events  are  crustal,  located  at  depths 
less  than  30  km.  Initially,  waveforms  were  analyzed  if:  (1)  the  event  was  within  4°  of  the 
network  so  that  P  and  5  waves  were  well  defined,  (2)  at  least  four  phases  were  recorded  at 
least  by  three  stations  for  the  location,  (3)  the  signal-to-noise  ratio  of  main  arrival  to  pre¬ 
event  noise  was  greater  than  3,  and  (4)  there  were  no  telemetry  gaps  in  recording.  Data  for 
1992  were  picked  and  associated  by  the  Joint  Seismic  Program  (JSP)  Center  (Harvey  et  al. 
1994).  For  1993  and  1994,  the  arrivals  were  picked  visually  using  JSP  software  and 
events  were  located  using  locSAT  (Bratt  &  Bache  1988).  A  total  of  249  P  and  304  S 
waveforms  were  used. 

Measurement  oft*.  Initial  observations  show  that  seismic  waves  propagate  efficiently 
through  platform  paths,  whereas  high  frequency  signals  become  weaker  for  paths  that 
cross  mountain  region  (Figure  2).  Hence,  events  used  here  were  classified  into  two 
groups:  those  with  ray  paths  that  sample  the  Russian  Platform  (plusses  in  Figure  1)  and 
those  with  paths  that  sample  the  Greater  Caucasus  mountains  (solid  circles  in  Figure  1). 
Frequency-independent  attenuation  parameters  are  measured  by  fitting  body  wave  spectra 
H(f)  (e.g.  Lees  &  Lindley  1994;  Lindley  &  Archuleta  1992;  Boatwright  1978)  to 

(» 
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where  H{f)  is  the  Fourier  displacement  amplitude,/is  the  frequency,  1(f)  is  the  known 
instrument  response  to  ground  displacement,  flo  is  the  frequency-independent  term,/c  is 

the  source  comer  frequency,  /*  is  the  frequency-independent  attenuation  operator.  The  P 
and  5  arrivals  are  selected  in  a  2-5  s  long  window  that  contains  the  dominant  energy,  and 
their  Fourier  amplitude  spectra  are  calculated  by  a  multitaper  technique  (Park  et  al.  1987).  S 
waves  are  rotated  to  the  transverse  orientation.  For  each  record,  a  time  window  is  chosen  to 
include  all  direct  arrivals  and  to  encompass  most  of  the  initial  arrival  energies.  Window 
length  typically  varies  from  2  to  4  s  for  P  waves  and  from  3  to  5  s  for  S  waves,  depending 

on  epicentral  distance.  Tests  showed  that  variations  of  t  with  window  length  are  small 
and  have  no  significant  affect  on  Q  measurements. 

The  log-amplitude  spectra  of  P  and  S  waves  are  fit  in  a  least  squares  sense  to  log[H(f)] 
to  determine  /*  ,fc  and  CIq  (Figure  3).  The  t  values  include  several  large  outliers.  In 
order  to  ensure  robustness  in  inversion,  we  eliminate  data  if  the  t  estimates  are  more  than 
3-0  from  the  mean  value  for  each  station.  We  assume  t  to  be  independent  of  frequency, 
probably  reasonable  over  the  relatively  small  (2-15  Hz)  frequency  of  most  of  our 
measurements.  For  both  P  and  S  records,  noise  levels  are  determined  from  the  spectrum 
of  a  window  placed  just  prior  to  main  arrivals.  These  noise  measurements  are  used  to 
determine  the  usable  frequency  band  for  each  signal,  limited  to  frequencies  where  the 
signal-to-noise  ratio  is  obviously  greater  than  1.  As  well,  frequencies  below  the  instrument 
natural  frequency,  0.2  Hz,  are  not  used.  Typical  frequency  ranges  are  for  P  waves  is  1-15 
Hz,  and  1-12  Hz  for  S  waves. 

In  the  initial  f  measurements,  the  fc  for  a  single  event  varies  widely  between  stations, 
perhaps  because  t‘  estimates  trades  off  with  fc-  In-order  to  understand  any  possible 
systematic  correlation  between  these  parameters,  t  is  estimated  twice:  once  by 
constraining/c  to  a  single  value  for  all  records  of  the  same  event  at  different  stations,  and 


once  by  allowing/^  to  vary  between  stations.  Results  are  similar  (Table  1),  indicating  that 
the  trade-offs  do  not  significantly  affect  our  results,  a  result  also  found  by  Lees  &  Lindley 
(1994).  In  all  subsequent  discussions,  we  present  Q  measurements  from  unconstrained /c 
case  which  are  based  on  more  observations. 


1.4  Path-averaged  Attenuation 


Initially,  we  separate  the  t  measurements  into  two  paths:  platform  paths  and  mountain 
paths.  For  each  subset  of  t*  measurements  we  estimate  a  path-averaged  attenuation 
operator,  Q  (e.g.  Hough  et  al.  1988a)  by 


*  X  \  * 

t  (r)  =  to 


(2) 


where  is  the  near-receiver  contribution  term,  r  is  the  distance.  Here,  V  is  the  path- 
averaged  velocity  of  relevant  phases;  6.0  and  3.5  km/s  for  crustal  P  and  5,  respectively 
(Balavadze  &  Tvaltvadze  1958).  The_mantle  Pn  velocity  is  assumed  to  be  8.0  km/s,  and 
4.6  for  Sn-  For  this  simple  model,  Q  and  r^are  determined  from  linear  regression  of  the 
i’  (r)  measurements  (Table  1).  We  find  that  Q  for  mountain_paths  is  two  to  three  times 
lower  than  platform  paths,  and  Q  for  P  and  5  waves  ( )  are  roughly  equal  (Figure 
4,  left). 

At  the  five  stations,  (q  varies  from  0.012  to  0.046  s  for  P  waves  and  0.016  to  0.051  s 
for  S  waves,  values  which  are  small  but  non  zero.  We^ interpret  this  parameter  as  the 
attenuation  caused  by  lossy  near-surface  rocks.  The  tg  values  are  systematically  larger  for 
stations  located  on  poorly  consolidated  sediments  of  the  Russian  Platform  (NAG,  MIC  and 
KUB)  than  for  mountain  stations  (KIV  and  GUM)  (Table  1  and  Figure  4,  right),  at  least 
for  S  waves. 


1.5  Inversion 

In  order  to  better  understand  the  variations  of  observed  f*  estimates,  we  modified  the 
velocity  inversion  technique  of  Roecker  (  1993)  to  invert  for  Q  ' .  The  t  estimates  are 
related  to  Q''  by 

/•=  f  (3) 

where  the  integral  is  taken  along  the  ray  path  of  seismic  wave.  Because  t  is  linearly  related 
to  Q~' ,  the  inversion  and  error  assessment  here  become  much  simpler  than  for  velocity 
inversion.  In  practice  we  discretize  the  integral  and  assume  that  both  Q  and  V  are  constant 
within  small  blocks. 

Parameterization  for  Q~' .  Rather  than  present  results  for  a  single  inversion  with  large 
numbers  of  free  parameters,  we  invert  t‘  measurements  for  Q~'  in  a  sequence  of 
increasingly  complex  parameterizations  to  test  the  significance  of  features  of  interest  (Table 
2).  In  all  inversions,  the  attenuation  grids  are  oriented  N20°E  so  that  one  coordinate  is 
parallel  to  the  strike  of  the  mountain.  The  medium  is  parameterized  by  blocks  with  constant 
(2"'  and  constant  velocity,  but  arbitrary  shape,  using  the  reparameterization  scheme  of 
Abers  &  Roecker  (1991).  Ray  paths  are  calculated  for  each  block  by  an  approximate  ray 
tracing  technique  similar  to  that  of  Thurber  &  Ellsworth  (1980).  In  some  inversions  vertical 
variations  in  Q''  are  included  by  parameterizing  velocity  into  three  layers:  a  near-surface 


layer  (200  m  thick),  crust  (45  km  thick  (Alexidze  et  al.  1993))  and  upper  mantle  (below  45 
km  depth).  In  others,  the  Earth  is  assumed  to  consist  of  a  single  layer  below  the  200  m 
thick  near-surface  layer.  The  near-surface  layer  is  included  to  mimic  site-effects  of  soft 
sediments  and  is  parameterized  such  that  a  separate  block  lies  beneath  each  station.  For 
both  P  and  S,  blocks  that  hit  by  at  least  four  rays  are  kept  for  inversion. 

Inversion  and  Constraints.  For  each  phase  {P  and  5),  Q  '  is  estimated  using  a  linear 
inverse  formalism  of  Tarantola  &  Vallette  (1982).  This  formalism  takes  into  account  a 
priori  uncertainties  of  data,  model  and  theory,  and  allows  a  convenient  and  realistic  way  to 
estimate  Q~'  along  with  its  a  posteriori  uncertainty  (formal  error)  and  resolution.  The 
solution  has  the  standard  form 

m  =  m„  -h [G\C,,  +  C„)-' G  +  C,]-' G\C,,  +  CJ"' -  Gm„) 

where  m  is  the  estimate  of  model  parameters,  mo  is  an  a  priori  estimate  of  m,  G  is  the 
matrix  of  partial  derivatives  of  t*  with  respect  to  model  parameters,  Cdd  is  the  a  priori 
covariance  matrix  of  the  data  vector  d,  Cmm  is  the  a  priori  covariance  matrix  of  m,  Ctt  is 
the  a  priori  covariance  matrix  of  the  theoretical  uncertainty.  Assuming  uncertainties  of  data, 
model  and  theoretical  errors  are  independent  and  uncorrelated,  the  corresponding 
covariance  matrices  Ccid>  Cmm  ^d  are  constructed  as  diagonal  with  individual  elements 
being  a] ,  cr^  and  erf ,  respectively. 

The  term  Om  and  mo  reflects  a  tie  between  the  resultant  model  and  a  priori  estimate  of 
the  model.  We  fix  mo  of  0.02  (0=50)  for  200  m  thick  near-receiver  layer  and  0.001 
(0=  1 000)  for  lower  layers.  The  a  priori  uncertainties  (Om)  are  set  equal  to  mo  in  each  ^ 
block.  The  a  priori  uncertainties  of  data  (oy)  are  set  to  the  standard  errors  calculated  for  t 
estimates,  and  are  typically  in  the  range  of  0.002  -  0.007  s  for  both  P  and  S  waves.  Also, 
we  fix  a/=0.01  s,  the  same  order  of  magnitude  as  the  r*  estimates,  as  a  crude  estimate  of 
uncertainty  due  to  incompleteness  of  the  parameterizadon. 

Source  depth  effects.  For  events  that  are  well  outside  the  network,  source  depth  is 
poorly  constrained.  When  coverage  is  good,  the  best  located  earthquakes  in  this  region  are 
no  deeper  than  -15  km  (Triep  et  al.  1995),  so  it  seems  reasonable  to  assume  that  all 
Caucasus  events  lie  in  the  upper  crust.  (Some  hypocenters  used  in  the  present  study  are 
placed  by  location  programs  at  -100  km  depth,  where  continental  lithosphere  is  believed  to 
have  little  strength.  These  events  are  several  hundred  km  outside  the  network  and  the 
depths  are  probably  erroneous.)  Although  source  depth  is  weakly  constrained,  for  events 
outside  of  the  network,  the  actual  depth  has  little  effect  on  0  measurements.  A  series  of 
inversions  (not  shown)  were  conducted  in  which  all  event  depths  are  constrained  1,  5,  and 
20  km.  Residual  variances  do  not  show  any  significant  change  between  these  inversions 
{P  <3%;  S  <1%),  and  0"'  values  are  similar  within  their  formal  uncertainties.  Hence, 
reasonable  depth  uncertainties  will  not  grossly  affect  0”’  estimates. 

1.6  Results 

The  increasingly  complex  sequence  of  inversions  (Table  2)  show  trade-off  between 
residual  variance  and  the  number  of  free  parameters  (Figure  5).  The  variance  decrease  is 
less  than  1%  between  final  inversion  shown  (parameterization  J,  Figure  6)  and  most 
complex  parameterization  attempted  (parameterization  K).  The  data  variance  in  J  is  0.609  x 
10-3  s2  for  5  and  0.341  x  lO’^  s^  for  P,  a  -60%  reduction  relative  to  the  homogeneous 
half-space  parameterization  (A).  For  all  inversions,  the  a  posteriori  error  incre^es  as  size 
of  block  decreases.  Uncertainties  associated  with  0;'  are  2-4x  larger  than  0;' .  Inversion 
results  for  0”'  along  with  a  posteriori  errors  are  shown  in  map  views  and  cross  sections 


for  two  different  parameterizations  H  and  J  (Figure  6  and  7).  For  crustal  blocks,  resolution 
diagonals  are  >0.9  ( Q~' )  or  >0.84  ( Q~' ).  For  mantle  blocks  ray  coverage  is  sparser,  and 
resolution  diagonals  vary  from  0.48  to  0.99  (Qj' )  or  0.36-0.96  (G^' )• 

From  these  tests,  two  aspects  of  parameterizations  emerge  as  being  responsible  for 
most  of  the  variance  reduction:  allowing  different  G~'  values  beneath  the  mountains  from 
those  beneath  the  platform  (parameterizations  C,  D,  F,  H-K),  and  including  the  near¬ 
surface  "site  effect"  terms  {E-K).  Inversions  that  allow  attenuation  to  vary  across  the  strike 
of  the  mountains  (e.g.  parameterization  F)  shows  that  signals  are  roughly  2-3  times  more 
attenuated  beneath  the  mountains  than  beneath  the  platform  region  (for  mountain  paths  G”' 
and  G7'  are  1.39  ±0.09  x  10'^  and  1.36  ±0.06  x  lO'^,  and  for  platform  paths  they  are 
0.58  d±).07  X  10-3  and  0.38  ±0.04  x  lO'^,  respectively),  and  that  Q~'  is  roughly  equal  to 
Q~'  in  each  region.  The  corresponding  residual  variances  of  t  are  reduced  by  20-30% 
from  parameterizations  that  require  lateral  homogeneity  (compare  inversions  C,  D,  and  F 
with  A,  B,  and  E  respectively).  For  both  P  and  5,  inversions  that  account  for  both  lateral 
and  vertical  variations  in  Q;  (parameterization  H  and  J)  show  that  attenuation  contrasts 
between  mountain  and  platform  paths  are  similar  to  parameterization  F,  similar  for  both 
crust  and  upper  mantle  (Figure  6a  and  b  (middle  and  bottom)  and  7).  Q  '  values  for  near- 
surface  station  blocks  vary  between  0.002  and  0.08,  with  the  highest  values  seen  at 
stations  on  soft  and  young  sediments  north  of  the  network  (Figure  6a  and  b,  top). 

The  total  variance  reduction  for  inversion  J  is  less  than  5%  relative  to  inversion  F,  the 
simplest  to  allow  both  near-station  effects  and  lateral  variation  between  mountains  and 
platform.  By  comparison,  variance  reduces  by  more  than  50%  between  inversions  F  and  A 
(homogeneous  half-space).  The  variance  decrease  is  small  between  inversions  that  allow 
for  separate  upper  mantle  structures  and  those  that  do  not,  6%  (compare  Model  G  and  H), 
implying  that  the  mantle  structures  are  not  well  resolved.  F-tests  show  that  a  lateral 
boundary  between  mountains  and  shield  is  required  by  the  data,  but  further  complexity  is 
not. 

Resolution.  A  synthetic  data  set  of  t‘  values  is  inverted  for  G  '  •  The  artificial  data 
are  calculated  for  a  model  with  a  single  lateral  step  (Figure  8,  top  left),  using  all  ray  paths 
and  uncertainties  in  the  real  data  set.  When  inverted  with  parameterization  {H),  the  correct 
G  '  is  recovered  to  within  10%  in  the  crust  and  14%  in  the  mantle  (Figure  8,  top  right). 
When  inverted  with  a  more  complex  parameterization  (7),  results  depart  from  the  correct 
values  by  up  to  18%  for  crustal  blocks,  although  the  overall  pattern  is  recovered  (Figure  8, 
bottom).  Mantle  blocks  vary  more,  up  to  54%.  The  variance  reduction  gained  by  adding 
mantle  blocks  is  less  than  10%,  because  most  of  the  ray  paths  are  unaffected  by  the  upper 
mantle  structures.  Resolution  diagonals  are  30-40%  lower  for  mantle  blocks  relative  to 
simpler  inversions.  For  all  blocks,  the  deviations  between  true  models  and  that  determined 
by  the  inversion  are  less  than  l-o  formal  errors.  This  suggests  that  the  formal  error 
estimates  may  overestimate  true  uncertainties  by  several  times,  although  several  sources  of 
error,  such  as  error  associated  with  incomplete  parameterization,  are  not  modeled  in  this 
test. 

Boundary  Between  High-Q  and  Low-Q  regions.  The  boundary  between  high-G  snd 
low-(2  regions  is  sharp  and  closely  follows  the  northern  foothills  of  the  mountains.  The 
location  of  this  boundary  is  evaluated  by  performing  inversions  repeatedly  for 
parameterizations  with  a  single  lateral  Q  step,  similar  to  parameterization  F,  but  moving  the 
location  of  the  boundary  north  and  south  (Figure  9b).  For  each  inversion,  the  residual 
variance  is  calculated.  Based  on  this  variance  reduction,  for  both  P  and  S  the  location  of  the 
boundary  is  constrained  at  43.67°E  longitude  to  lie  at  43.265°N  latitude,  along  the  lower 
foothills  of  the  northern  slope  of  the  Greater  Caucasus.  The  95%  confidence  limits  (F  test) 


on  N-S  boundary  location  are  80  km  wide  for  S  and  130  km  wide  for  P  waves  (Figure  9a). 
Hence  the  boundary  is  beneath  the  northern  foothills,  where  surface  geology  shows  a 
simple  monoclinal  fold. 

Depth  sampling.  Although  the  inversion  results  show  similar  Q  for  the  crust  and 
upper  mantle  beneath  the  mountains  and  beneath  the  platform,  it  is  unclear  whether  the 
present  data  can  resolve  upper  mantle  attenuation  independently  of  that  in  the  crust.  Ray 
tracing  shows  that  only  a  few  ray  paths  (~12%)  sample  the  upper  mantle  (although  source 
depths  are  uncertain).  The  6%  variance  reduction  between  inversions  that  include  a  separate 
crust  and  mantle  (inversions  H,  J)  and  those  that  do  not  (inversions  F,  /)  is  statistically 
insignificant  above  a  90%  confidence  level,  so  Q  variations  with  depth  are  not  constrained 
by  the  present  data  set.  The  point  is  further  illustrated  by  examining  the  windows  chosen 
for  initial  r*  measurements.  Examining  the  waveforms  for  paths  between  1.8  and  3.5 
distance  (Pn-Pg  cross-over  distance  is  -1.8°)  shows  that  the  selected  time  windows  include 
both  crustal  (fI  Pb)  and  mantle  (Pn)  phases  (Figure  10).  Spectral  measurements  will  be 
dominated  by  me  larger  signal,  which  is  typically  Pg.  Hence,  it  is  reasonable  to  infer  that 
the  variations  in  Q~'  are  best  considered  vertical  averages  of  me  entire  crust. 

1.7  Discussion 

Attenuation  values.  Measurements  made  here  are  in  the  2-15  Hz  range,  and  can  be 
compared  to  values  determined  in  omer  regions  at  these  frequencies.  The  Q  values  for  the 
Russian  shield  (Table  1)  are  similar  to  mose  made  in  other  stable  cratons,  while  values 
from  the  Caucasus  are  typical  of  tectonically  active  regions.  For  example,  several  studies 
show  that  Qs  for  eastern  North  America  is  near  2000  at  1-10  Hz  frequencies  (Atkinson  & 
Mereu  1992;  Boatwright  1994;  Frankel  et  al.  1990),  similar  to  the  Russian  shield  values 
found  here.  Measurements  from  California  show  Qp  close  to  Qs,  near  800  at  Loma  Prieta 
(Fletcher  &  Boatwright  1991)  or  near  1000  at  Anza  (Hough  et  al.  1988a).  The  California 
values  are  similar  to  our  estimates  for  attenuation  beneath  the  Caucasus.  A  preliminary 
estimate  for  Q  in  the  Kopet  Dagh,  near  the  Turkmenistan-lran  border,  shows  Qp~  Qs  ~ 
700.  somewhat  lower  than  values  found  in  the  Caucasus  (Abers  &  Sarker  1996). 
Somewhat  lower  coda  Q  values  for  the  Caucasus.  Qcoda  =  330  (Rautian  et  al.  1979),  are 
made  from  1  Hz  coda  measurements  and  may  reflect  differences  in  frequency  and  wave 
type. 


I.7.i  Attenuation  Mechanism.  In  general,  the  causes  of  apparent  attenuation  are 
complex  and  include  scattering  by  crustal  heterogeneity,  presence  of  cracks  that  are  partially 
or  fully  saturated  with  fluids,  and  intrinsic  dissipation  (e.g.  Mitchell  1995;  Aki  &  Chouet 
1975;  Hough  &  Anderson  1988b;  Frankel  1982).  In  many  situations  it  is  unclear  that 
scattering  and  intrinsic  mechanisms  can  be  separated.  However,  some  diagnostics  favor 
intrinsic  mechanisms  for  our  observations.  Overall,  any  feasible  explanation  of  our  data 
should  account  for  several  observations:  (I)  2”'  and  Q^  are  roughly  equal  (±11%) 
beneath  the  Greater  Caucasus  and  beneath  the  Russian  Platform;  (2)  for  both  P  and  S,  Q 
for  mountain  paths  is  2-3  times  larger  than  that  of  platform  paths;  and  (3)  a  sharp  boundary 
between  high-2  ^^d  low-2  follows  the  northern  foothills  of  the  mountains. 

Scattering.  When  a  body  wave  propagates  through  a  heterogeneous  medium  some 
high-frequency  energy  is  lost  to  the  wavefield  coda  by  scattering  (e.g.  Wennerberg  & 
Frankel  1989).  Both  scattering  and  intrinsic  absorption  contribute  to  observed  Q  of  body 
waves.  Energy  conservation  arguments  and  numerical  simulations  show  that  in  a  medium 
dominated  by  scattering,  coda  will  be  enriched  in  high  frequencies  relative  to  body  waves 


(Richards  &  Menke  1983;  Frankel  &  Wennerberg  1987).  By  comparison,  intrinsic 
attenuation  represents  a  loss  of  energy  at  high  frequencies  throughout  the  seismogram. 

To  test  for  scattering,  we  compare  the  Fourier  amplitude  spectra  of  coda  signals  to 
those  of  the  main  body  wave  pulses  for  a  small  subset  of  events  used  in  this  study  (Figure 
1 1).  Early  coda  spectra  are  measured  in  a  time  window  that  is  10-15  s  delayed  from  the 
corresponding  body  wave  arrival,  and  late  coda  spectra  are  measured  at  twice  the  5-wave 
travel  time.  The  sharp  spectral  fall-offs  seen  for  body  waves  sampling  mountain  paths  are 
also  seen  the  corresponding  coda  spectra,  with  similar  spectral  slopes.  Such  sharp  high- 
frequency  decays  are  not  seen  for  shield  paths.  Hence  it  appears  that  both  direct  and  coda 
signals  propagating  through  mountains  are  attenuated  similarly  at  high  frequencies,  a 
characteristic  of  intrinsic  attenuation  but  not  obviously  a  characteristic  of  scattering.  It  is 
less  clear  that  Q  for  shield  paths  are  dominated  by  intrinsic  attenuation,  as  spectral  fall-off  is 
small  below  -8  Hz  for  these  paths,  and  fall-off  in  some  cases  is  sharper  for  body  waves 
than  for  coda  (Figure  1 1).  Hence,  body-wave  Q  values  (~2000)  seen  for  the  Russian  shield 
conceivably  could  be  dominated  by  scattering,  but  the  low  Q  seen  for  mountain  paths  is 
more  consistent  with  intrinsic  attenuation. 

Focusing.  Q  is  determined  here  from  spectral  fall-off,  so  the  only  focusing  effects  of 
interest  are  those  that  produce  variation  in  the  2-15  Hz  frequency  band.  This  requires 
structures  that  are  of  a  size  comparable  to  the  relevant  wavelengths  (0.2  -  10  km). 

However,  similar  attenuation  patterns  are  seen  at  all  stations,  over  distances  of  50-100  km, 
and  sources  within  each  tectonic  province  are  distributed  over  several  hundred  km.  The 
same  focusing  structure  cannot  be  responsible  for  producing  frequency-dependent  effects 
along  all  of  these  paths.  Some  large-scale  amplitude  anomalies  may  be  attributed  to 
focusing  effects  of  Moho  structure  (Cormier  &  Anderson  1996). 

Cracks.  At  appropriate  temperatures  and  degrees  of  saturation,  the  presence  of  cracks 
can  lead  to  similar  Q,  and  Q~'  (e.g.  Mavko  &  Nur  1979),  agreeing  with  our  observation 
(1).  These  effects  may  be  responsible  for  Q''  patterns  seen  in  some  fault  zones  (Lees  & 
Lindley  1994).  However,  it  is  unlikely  that  cracks  will  remain  open  at  depths  greater  than 
6-8  km  (equivalent  to  -0.20  GPa)  (e.g.  Jackson,  Peterson  &  Fitz  Gerald  1992;  Christensen 
&  Wepfer  1989).  Most  of  the  ray  paths  studied  here  are  >100  km  long  and  spend 
substantial  time  in  the  deeper  parts  of  the  crust  or  upper  mantle,  where  open  unsaturated 
cracks  are  unlikely.  Cracks  may  contribute  to  attenuation  within  the  near-surface  layers 
(<200  m  depth). 

Anelastic  Absorption.  Most  mechanisms  of  intrinsic  attenuation  depend  upon 
temperature  (e.g.  Karato  &  Spetzler  1990;  Anderson  1989),  and  high  temperatures  or 
partial  melts  have  been  considered  as  primary  attenuation  mechanism  beneath  several 
mountain  belts  including  the  Caucasus  (e.g.  Rodgers  et  al.  1997).  Of  the  few  laboratory 
studies  that  have  been  made  at  high  temperatures  and  seismic  frequencies,  high-pressure 
experiments  on  dunites  (Jackson  et  al.  1992)  and  room  pressure  experiments  on  basic 
rocks  (Kampfmann  &  Berckhemer  1985)  both  show  that  Qs  and  Qp  drop  by  a  factor  of  5- 
10  as  temperature  increases  close  to  the  solidus.  Hence,  large  Q  variations  are  viable 
without  the  presence  of  partial  melts.  At  temperatures  well  below  the  solidus  (e.g.  <1000 
C  in  dunites)  attenuation  appears  to  operate  by  a  somewhat  different  mechanism  than  at 
higher  temperatures,  and  shows  an  approximately  linear  (non- Arrhenius)  relation  between 
logi  Q~' )  and  temperature  (Kampfmann  &  Bercldiemer  1985;  Jackson  et  al.  1992). 
Although  individual  studies  show  large  differences  in  parameter  values  (see  next  section), 
in  each  case  a  factor  of  three  attenuation  variation  as  seen  here  can  be  easily  explained  by 
temperature  variations  well  below  the  solidus. 


Our  observation  that  Qp/Qs  ~  1  differs  from  most  deep-mande  observations,  which 
show  Qp  /Qs~  2  (e.g.  Karato  &  Spetzler  1990).  Observations  of  Qp/Qs  -  2  require  that 
attenuation  in  shear  {Q^r^)  be  much  larger  than  attenuation  of  bulk  modulus  {Qk'^), 
consistent  with  several  dissipative  mechanisms  (defect-based  or  grain-boundary  based) 
likely  to  occur  at  high  temperature.  By  contrast,  Qp/Qs~  ^  requires  that  Qfr^  and  Qk'^  be 
similar  in  size.  The  most  likely  mechanism  to  produce  large  Qk'^  is  thermoelasticity,  which 
is  efficient  in  polycrystalline  materials  with  grains  that  are  anisotropic  and  have  strongly 
varying  thermal  properties  (Anderson  1989).  It  is  easy  to  imagine  that  such  conditions 
exist  within  metamorphic  rocks  common  to  continent^  collision  zones,  although  little 
experimental  evidence  exists  that  quantifies  Qk'^  and  in  such  rocks.  As  experimental 
measurements  become  available,  the  joint  observation  of  Qp  and  Qs  then  may  have  value  in 
constraining  lower  crustal  composition. 

L7.ii  Temperature  beneath  the  Greater  Caucasus.  If  the  observed  Q 
variations  between  shield  and  Caucasus  mountains  are  entirely  due  to  temperature,  then 
these  observations  can  constrain  the  temperature  difference.  At  temperatures  below  900°  C, 
experiments  on  mafic  and  ultramafic  rocks  suggest  that  temperature  (7)  is  related  to 
attenuation  by 

r  =  1100-150/^“/og,o((2,)±50°C  (5) 

(Kampfmann  &  Berckhemer  1985).  Using  this  formula,  our  measurements  made  at 
frequencies  of  2-15  Hz  (Table  1)  indicate  temperatures  of  520-580°  C  for  shield  paths  and 
590-650°  C  for  mountain  paths,  for  a  temperature  contrast  of  only  -70°  C  at  any  frequency 
(based  on  5  wave  1000/2  values  with  unconstrained  comer  frequencies.  Table  1).  The 
shield  temperatures  are  reasonable  for  values  near  Moho  depths  in  old  continents  (e.g. 
Sclater,  Jaupart  &  Galson  1981),  consistent  with  our  inference  that  the  waves  sampled 
mostly  travel  in  the  lower  crust  and  uppermost  mantle.  The  small  temperature  increase  is 
insufficient  to  cause  melting  in  mafic  or  ultramafic  rocks. 

Two  factors  temper  our  confidence  in  these  numbers.  First,  the  Q  values  for  shield 
paths  may  be  lowered  by  scattering  (see  above)  and  may  represent  an  upper  bound  on 
temperature  there.  Second,  the  experiments  made  by  Kampfmann  and  Berckhemer  were 
made  at  room  pressure,  not  at  lower  crustal  conditions.  In  experiments  at  high  pressures, 
Jackson  et  al.  ( 1992)  show  a  much  larger  factor  of  ten  increase  in  2”'  over  the  temperature 
interval  0°  -  1000°  C  for  dunite.  Although  they  show  that  Qs  will  equal  775  at  temperatures 
below  600°  C,  a  factor  of  2.6  variation  in  Qs  (Table  1)  would  translate  into  a  -400°  C 
temperature  contrast  between  shield  and  mountains,  about  5  times  larger  than  predicted  by 
Kampfmann  and  Berkhemer's  experiments.  Since  the  experiments  of  Jackson  et  al.  (1992) 
were  made  in  dunite  and  concentrated  on  near-solidus  mechanisms  their  applicability  to  our 
setting  is  uncertain. 

In  conclusion,  moderate  temperature  increases  beneath  the  Greater  Caucasus,  sufficient 
to  produce  regional  metamorphic  effects,  may  explain  the  attenuation  contrast.  Although 
any  direct  association  with  temperature  is  not  yet  well  constrained,  small  temperature 
variations  are  most  consistent  with  our  observations.  The  presence  of  Plio-Quatemary 
volcanism  in  parts  of  the  Greater  Caucasus  shows  that  some  melting  must  be  occurring  at 
depth,  but  it  need  not  be  widespread  nor  be  the  primary  cause  of  the  attenuation  variations. 
Although  most  previous  studies  do  not  give  quantitative  Q  values,  it  seems  likely  that 
previous  observations  of  strong  S„  attenuation  beneath  the  Caucasus  (e.g.  Kadinsky-Cade 
et  al.  1981;  Rodgers  et  al.  1997)  can  be  explained  by  elevated  temperatures  and  not  partial 
melt. 


I.T.iii  Geological  Implications.  Our  results  show  high  attenuation  beneath  high 
topography,  with  a  sharp  boundary  that  closely  follows  the  northern  foothills  of  the 
mountain  ranges.  Geologically,  the  boundary  is  represented  by  the  active  deformation  front 
lies  south  side  of  the  mountain  ranges  (e.g.  Philip  et  al.  1989;  Triep  et  al.  1995; 

Zonenshain  &  Le  Pichon  1986).  Surface  geology  also  shows  that  the  thrust  sheets  that 
characterize  deformation  on  the  southern  slope  are  transported  further  south  over  adjacent 
massif  (e.g.  Dotduyev  1986;  Philip  et  al.  1989).  The  inferred  boundary  at  depth,  thus,  is 
north  of  the  structural  boundary  at  the  surface,  and  underlies  a  region  where  surface 
geology  shows  only  a  simple  monocline.  In  other  words,  high  attenuation  ((2“')  is 
observed  not  just  where  faults  are  seen  at  the  surface,  but  at  depth  where  they  project  down 
dip. 

If  our  inferences  for  the  mechanism  of  Q~'  are  correct,  then  some  process  is  heating  the 
mountain  range  throughout  its  extent.  One  speculative  possibility  is  that  convective 
instability  within  the  mantle  is  elevating  temperatures  as  a  consequence  of  lithospheric 
thickening  (e.g.  Houseman,  McKenzie  &  Molnar  1981).  Such  a  mechanism  would  help 
explain  the  correspondence  between  inferred  high  temperatures  and  topography,  unlike 
scenarios  that  appeal  to  high  temperatures  remnant  from  pre-Miocene  subduction.  Another 
possibility  is  that  shear  heating  is  significant  at  the  downdip  end  of  the  active  Greater 
Caucasus  thrust  fault  systems  (e.g.  Figure  1),  either  through  brittle  faulting  or  through 
viscous  dissipation  in  ductile  shear  zones.  If  the  low  P„  velocities  observed  by  Hearn  and 
.Ni  ( 1 994)  for  the  Caucasus  as  a  whole  coincide  with  the  attenuation  anomalies  seen  here, 
then  that  would  suggest  that  the  region  of  high  temperatures  extends  too  deep  to  be 
attributable  to  shear  heating  on  faults.  More  work  on  the  depth  extent  of  structures  needs  to 
be  done  to  test  these  possibilities. 

High  attenuation  beneath  the  Greater  Caucasus  would  cause  densities  to  be 
anomalously  low.  As  a  consequence,  some  component  of  the  high  topography  may  be 
supported  by  Pratt-type  mechanisms  rather  than  by  Moho  deflections  as  is  often  assumed 
(Philip  et  al.  1989;  Ruppel  &  McNutt  1990).  For  example,  a  temperature  increase  of  70  - 
400°  C  throughout  a  150  km  thick  lithosphere  would  generate  0.3  -  1.8  km  of  topography, 
for  a  coefficient  of  thermal  expansion  of  3  x  lO’^  ’  K‘‘.  One  consequence  is  that  the 
Bouguer  gravity  low  over  the  Caucasus  may  not  be  a  very  good  proxy  for  Moho  depth 
changes.  Unfortunately  we  know  of  no  modem,  direct  measurements  of  cmstal  thickness 
beneath  the  Caucasus  to  test  this  speculation. 

L7.iv  Implications  for  Regional  Monitoring.  No  matter  what  the  cause  of 
Q~' ,  the  variations  in  signal  behavior  between  mountain  and  shield  paths  has  direct 
implications  for  use  of  body  wave  amplitudes  in  source  discrimination  and  yield  estimation. 
Similar  variations  in  Lg  and  Sfj  amplitudes  are  seen  for  longer  paths  in  the  region  (Abers, 
Kim  &  Lemer-Lam  1995;  Rodgers  et  al.  1997).  The  sharp  boundary  between  high- 2"' 
and  low-  Q~'  regions  suggests  that  the  stations  placed  on  the  mountains  suffer  more 
attenuative  loss  of  signal  than  those  on  the  nearby  shield.  More  importantly,  smooth 
attenuation  models  such  as  those  based  upon  long-period  surface  waves  probably  alias 
attenuation  changes  across  major  tectonic  boundaries.  Such  sampling  problems  are 
significant  for  paths  that  lie  within  a  few  hundred  km  of  the  boundary.  Propagation  across 
such  short  paths  could  play  an  important  role  in  verifying  the  Comprehensive  Test  Ban 
Treaty,  as  the  characteristics  of  signals  generated  by  small  events  come  under  increasing 
scrutiny. 
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11.  Long  Path  Data  Analysis. 

Summary.  As  part  of  this  project  we  extract  a  set  of  waveforms  suitable  for  long-path 
analysis  from  the  Caucasus  Network  (CNET)  data  set.  Because  the  primary  recording  was 
on  short  (~3  minute)  triggered  windows,  late  phases  were  missed  on  many  seismograms 
for  path  lengths  >5°  and  some  events  were  not  recorded  at  all.  A  secondary,  backup  set  of 
continuous  data  tapes  were  uncovered  and  decoded  to  provide  long  records  for  these 
events.  The  strategy  was  to  search  for  records  corresponding  to  all  Iranian  Plateau  events 
identified  in  the  PDE  catalog.  From  the  time  period  of  operation,  20  events  over  6  months 
of  recording  were  in  the  PDE  catalog  and  derivable  in  this  way  (Figure  12).  These  records 
demonstrate  several  aspects  of  the  Caucasus  network  recording  and  contemporaneous 
Geyokcha  array  recording.  In  particular,  (1)  the  CNET  triggering  algorithm  only  "missed 
1  of  20  events  that  had  visible  signals  at  >0.2  Hz,  due  to  a  month-long  triggered  data 
outage;  (2)  Lg  is  absent  and  Sn  is  only  present  at  frequencies  <  1  Hz  on  paths  >  1000  km  at 
both  CNET  and  Geyokcha;  (3)  detection  is  poor  for  Zagros  events  with  magnitudes  <4.5 
from  either  of  these  sites  (10  - 18*  range);  and  (4)  surface  waves  at  <  0.1  Hz  may  provide 
additional  useful  signal.  These  results  suggest  that  detection  of  small  events  is  difficult 
along  these  paths,  and  that  Sn/Pn  discriminants  here  need  to  be  carefully  designed  around 
the  optimal  frequencies  of  propagation. 

11. 1  Caucasus  Network 

The  Caucasus  network  was  established  in  1991  through  the  Joint  Seismic  Program,  as  a 
cooperative  project  between  Lamont-Doherty  and  the  OME-Obninsk  group  (Abers,  1994). 
The  network,  based  in  Kislovodsk  (KIV),  consists  of  5-6  digitally  telemetered  3- 
component  sensors  with  0.2  Hz  natural  periods,  digitized  at  60  sps  (Figure  1).  These 
stations  are  on  the  north  side  of  the  Greater  Caucasus  range,  and  straddle  the  transition 
from  the  Russian  Shield  to  the  Caucasus  collisional  zone.  The  northern  stations  lie  on  flat- 
lying,  poorly  consolidated  Neogene  sediments,  while  southern  stations  lie  on  tilted,  well- 
lithified  Mesozoic  sediments.  The  response  of  the  Kinemetrics  SV/SH  sensors  are  well- 
suited  for  analysis  of  local  and  regional  seismic  waves,  although  prior  to  mid- 1993 
triggering  windows  usually  did  not  include  S  or  Lg  for  events  farther  than  10-15  from  the 
network.  We  extracted  seismograms  from  a  continuously-recording  backup  system  in  order 
to  produce  a  waveform  data  set  suitable  for  these  longer  distances  (Figure  13).  Data  for 
1992  were  picked  and  associated  by  the  Joint  Seismic  Program  Center  (Harvey  et  al., 

1994)  and  subsequent  seismograms  are  distributed  through  IRIS.  A  total  of  645  events 
were  associated  for  1992  (292  in  the  PDE),  and  327  were  associated  with  the  PDE  catalog 
in  1993. 

11.2  Long-duration  records  for  Zagros  Events 

Continuous  recording  backup  tapes  were  generated  at  CNET  from  October  1993  to  May 
1994,  as  part  of  a  network  enhancement  prior  to  shutdown.  During  this  period  20  events 
occurred  in  southern  Iran  (in  the  Zagros  mountains)  that  were  listed  in  the  PDE  catalog. 
These  events  form  the  basis  for  this  study. 
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Long-path  records  were  retrieved  from  the  continuous  data  stream  for  these  20  events. 

Of  these  events,  5  triggered  the  system  and  so  the  P  waves  were  previously  detected 
(Figure  14).  The  trigger  algorithm  appeared  reasonably  robust,  as  only  one  PDE  events 
not  included  in  the  triggered  data  set  showed  any  visible  signal  in  the  continuous  data  (> 

0.2  Hz,  in  several  frequency  bands).  That  one  event  occurred  during  a  39  day  period  that 
the  triggered  stream  was  not  archived  (2/14/94  -  3/25/94).  Thus,  the  triggering  algorithm 
did  not  miss  any  visible  events. 

For  comparison,  these  same  20  events  were  compared  against  the  triggered  data  stream 
at  the  Geyokcha  array  in  Turkmenistan.  This  array  consisted  of  12  broadband  sensors 
deployed  over  a  ~3  km^  area,  with  supplemental  geophone  strings,  locates  in  the  Kopet 
Dagh  foothills  just  north  of  the  Iran  border  [Pavlis  et  al,  1994].  The  triggered  stream  only 
included  one  of  the  20  events,  with  comparable  signal-to-noise  characteristics  as  observed 
at  CNET.  Following  this  observation  we  suspect  that  the  observations  at  C^T  would 
apply  to  any  network  or  array  on  the  northern  margins  of  the  tectonically  active  collision 
belts  of  central  Eurasia,  although  quantitative  observations  have  been  made  only  at  this  one 
place. 

These  records  confirm  several  patterns  in  wave  propagation  identified  previously,  and 
quantify  them  further.  First,  the  Lg  wave  train  is  absent  on  paths  crossing  the  Caucasus 
(although  other  studies  show  it  present  at  the  same  stations  for  cratonic  paths).  Second,  Sn 
is  present,  but  only  at  frequencies  of  0.5  Hz  and  lower.  By  comparison,  previous  workers 
report  Sn  being  absent  [Kadinsky-Cade  et  al.,  1981;  Rogers  et  al,  1997].  We  suspect  that 
those  observations  are  a  frequency-dependent  property  of  wave  propagation  along  these 
paths,  and  at  least  for  the  older  studies  a  consequence  of  the  1  Hz  natural  frequency  of 
analog  Short-Period  sensors.  In  fact,  at  0.25  Hz  Sn  is  the  dominant  phase  on  the 
seismogram  and  Pn  is  only  marginally  visible  (Figure  13).  ^ 

A  final  observation  is  that  only  events  larger  than  mb  4.5  are  detectable  at  10-15  from 
CNET  (Figure  14).  This  observation  is  not  based  upon  a  great  number  of  points,  but  is 
consistent  with  the  catalog  analysis  for  1992  made  by  Harvey  et  al  [1994]  that  suggests 
CNET  under-detects  the  PDE  at  distances  greater  than  5°  from  the  network.  Again, 
comparison  with  records  at  the  Geyokcha  array  suggest  a  comparable  or  slightly  greater 
amount  of  attenuation  for  paths  to  the  Kopet  Dagh,  so  we  suspect  that  such  attenuation  is  a 
common  attribute  of  paths  that  traverse  the  Iranian  Plateau.  As  a  consequence  monitoring 
along  these  paths  alone  will  likely  be  insufficient  to  attain  CTBT  detection  goals. 

In  summary,  it  appears  that  events  of  magnitude  <  4.5  in  southern  and  eastern  Iran 
cannot  be  routinely  detected  near  the  Caucasus  network.  We  confirm  the  absence  of  Lg 
along  these  paths,  and  better  describe  the  behavior  of  Sn.  The  Sn  phase  group  does  exist 
but  only  at  frequencies  of  0.5  Hz  or  less.  Discriminants  based  on  Sn/Pn  amplitude  ratios 
may  be  possible  for  these  paths  but  they  would  require  careful  attention  paid  to  frequency 
of  measurement. 

III.  Amplitude  fluctuations  along  intermediate  paths 

Summary.  Amplitudes  of  seismic  waves  recorded  at  the  Caucasus  network,  along  the 
north  flanks  of  the  Greater  Caucasus,  are  measured  and  analyzed  for  attenuation 
characteristics.  From  one  year  of  observation,  96  events  between  1  and  10  from  the 
network  provide  stable  measures  of  RMS  Pn,  Sn,  Lg,  and  late  coda  amplitudes. 
Measurements  were  taken  from  seismograms  filtered  at  several  narrow  frequency  bands 
centered  from  0.5  to  8.0  Hz,  where  signal  levels  are  highest.  Our  results  confirm 
previously-inferred  spatial  variations  in  Sn  and  Lg  attenuation,  that  the  Greater  Caucasus 
marks  an  abrupt  boundary  between  the  high-2  Russian  Platform  and  a  region  of 
exceedingly  poor  Sn  and  Lg  propagation  within  the  collision  belt.  Paths  that  cross  large 
Quaternary  volcanic  provinces,  along  the  Greater  Caucasus,  seem  most  affected.  Amplitude 


ratios  show  the  largest  regional  differences  in  the  1-2  Hz  range  and  decrease  at  higher 
frequencies,  and  indicate  complicated  changes  in  the  mechanism  of  attenuation  between 
shield  and  tectonic  paths. 

For  a  given  path  the  RMS  amplitudes  of  the  Lg  group  and  late  coda  phases  predict 
magnitudes  as  well  as  can  be  expected  from  mb  uncertainties.  However  biases  of  1 .0-1.5 
magnitude  units  are  seen  for  RMS  amplitudes  along  different  paths,  at  frequencies  higher 
than  1  Hz.  Variation  is  much  reduced  by  taking  3-component  and  network  averages  of 
amplitude  measurements. 

111.1  Introduction 

Severe  changes  in  wave  behavior  are  observed  between  paths  that  cross  the  Caucasus 
and  other  Alpine-Himalayan  mountain  belts,  and  those  that  do  not  (Figure  15).  These 
changes  represent  substantial  attenuation  of  phases  often  used  to  discriminate  explosions 
from  earthquakes  and  to  estimate  event  size  and  type,  such  as  Lg  and  other  high-frequency 
surface-related  phases.  In  this  study  we  analyze  phase  amplitudes  for  events  within  10°  of 
the  Caucasus  Seismic  Network  (CNET).  Other  facilities,  such  as  the  Geyokcha, 
Turkmenistan  array  also  provide  digital  seismograms  that  give  complementary  information. 

111.2  RMS  Amplitudes  of  Regional  Phases 

We  have  measured  RMS  amplitudes  of  P-group,  5-group,  Lg  and  similar,  and  late  coda 
from  recordings  of  events  1-10°  from  CNet.  Because  the  CNet  time  series  contain 
occasional  spikes  and,  at  some  stations,  frequent  telemetry  gaps,  routine  and  automated 
spectral  methods  are  often  problematic.  More  robust  amplitude  estimates  are  made  from 
RMS  amplitudes  of  multiple-narrow-band  signals.  For  each  seismogram,  a  series  of 
narrow-band  filters  are  constructed  and  applied  centered  at  0.5,  1.0,  2.0, 4.0,  and  8.0  Hz. 
For  each  channel  and  frequency,  RMS  amplitudes  are  made  in  several  time  windows:  pre¬ 
event  noise,  Pn,  Sn,  Lg,  and  coda.  Pn  and  Sn  windows  extend  from  4  s  before  the 
predicted  phase  onsets  to  0.25  times  the  S-P  time  past  the  onset,  up  to  25  s  in  total  length. 
The  so-called  Lg  window  covers  the  period  in  group  velocity  between  3.6  and  2.6  km/s; 
the  range  extends  to  slower  than  normal  group  velocities  to  include  the  dominant  slow 
short-period  surface  phases  observed  for  shield  paths  (Figure  15).  The  coda  window  is 
centered  at  twice  the  5  travel  time  after  the  origin  lime,  and  is  40  s  in  duration.  Only 
segments  that  were  more  than  4  s  long  are  used,  and  only  those  where  dropouts  constituted 
less  than  10%  of  the  time  window.  These  requirements  allow  some  data  where  short 
trigger- windows  or  occasional  dropouts  might  be  present,  and  because  the  RMS 
measurements  are  integrated  quantities  the  measurements  are  still  valid.  All  seismograms  ■ 
are  visually  culled  for  other,  unusual  problems. 

An  RMS  combination  is  made  of  measurements  from  each  3-component  seismogram 
set,  for  as  many  components  as  were  usable.  Signal-to-noise  estimates  are  made  on  these 
combined  sets,  and  records  with  low  (<2)  signal  level  signals  are  eliminated.  Noise 
estimates  are  made  from  each  filtered  seismogram,  where  possible,  or  instead  taken  as  the 
maximum  observed  RMS  noise  level  for  a  particular  channel  and  frequency.  Amplitudes 
are  corrected  in  an  RMS  sense  for  incoherent  noise  contamination.  Finally,  amplitudes  are 
averaged  geometrically  across  the  network  for  each  event  (and  across  all  events  for  each 
station)  to  give  a  network-average  set  of  amplitude  ratio  estimates.  Only  averages  of  three 
more  stations  are  kept.  These  averages  were  considerably  less  variable  than  single  station 
estimates. 


III.3  Results 


In  map  view  (Figure  16)  large  variations  in  amplitude  ratios  are  seen  between  paths  that 
cross  the  crest  of  the  Greater  Caucasus  and  those  that  do  not,  and  indicate  extensive  Lg 
blockage  and  Sn  attenuation  throughout  the  actively  deforming  parts  of  the  Caucasus  belt. 
Such  behavior  was  seen  by  Kadinsky-Cade  et  al.  [1981],  although  they  did  not  have  access 
to  stations  north  of  the  Caucasus  and  were  limited  to  analog  records.  Using  some  of  the 
same  Caucasus  Network  data  Rodgers  et  al.  [1994, 1997]  saw  a  region  of  poor  S  and  Lg 
propagation  beneath  the  Caucasus.  Our  results  document  the  frequency  dependence  of  this 
phenomenon,  and  (similar  to  Rodgers  et  al.)  show  that  the  boundary  of  the  .attenuative 
region  is  sharp  and  lies  beneath  the  core  of  the  Greater  Caucasus.  Results  are  most 
consistent  for  the  Lg  window  although  the  Sn  window  behaves  similarly.  Differences  are 
most  pronounced  at  lower  frequencies  (Figure  17). 

Paths  from  events  in  the  Dagestan-Chechnya  region  (DA,  Figure  17, 18)  are  roughly 
parallel  to  the  strike  of  the  ranges  and  follow  the  margin  of  the  mountains  and  traverse  thick 
foredeep  sediments  north  of  the  Greater  Caucasus.  These  paths  are  also  characterized  by 
slow  (-2.9-3.0  km/s)  high-frequency  Rayleigh  waves  (Figure  15).  Large  amplitudes  are 
seen  most  frequencies,  particularly  from  0.5  to  2  Hz,  and  may  reflect  the  a  layered, 
relatively  undisrupted  nature  of  the  paths. 

Paths  that  cross  the  Greater  Caucasus  show  very  low-amplitude  S  and  Lg-related  phases 
relative  to  the  P  arrivals,  even  for  relatively  short  paths.  Many  events  are  roughly  200  km 
distant,  near  the  1991  Racha  earthquake  zone  (Triep  et  al.,  1995),  yet  show  amplitude 
ratios  3-10  times  lower  than  for  Dagestan  paths.  Waves  from  these  events  traverse  the 
active  deformation  front  of  the  Greater  Caucasus,  and  the  region  of  Plio-Quatemary 
volcanism  near  Mt.  Elbruz  at  the  SE  edge  of  the  network.  Paths  from  events  farther  south 
that  traverse  the  Lesser  Caucasus  show  very  small-amplitude  late  phases. 

Absolute  RMS  amplitudes  for  a  given  region  correlate  well  with  eveat  size  measures 
such  as  mb  (Figure  18)  although  the  scaling  parameters  vary  systematically  between 
regions.  The  mb  values  come  from  the  PDE  for  region  TU,  and  from  the  CNet  Information 
Product  for  closer  events  (so  that  some  coupling  with  RMS  Pn  is  expected).  In  general 
correlation  coefficients  are  0.8-0.9  for  each  population  between  1  and  4  Hz.  Events  from 
Turkey  show  RMS  amplitudes  of  late  phases  a  factor  of  10-100  low  at  2  Hz  compared  with 
shield  paths,  with  the  discrepancy  increasing  with  frequency.  Some  of  the  best  correlations 
(>0.9)  come  from  late  coda  amplitudes,  which  have  long  been  used  for  magnitude 
determination  in  the  Caucasus  (Rautian  et  al.,  1979).  The  variation  from  region  to  region 
corresponds  to  biases  of  1.0- 1.5  magnitude  units  in  the  ability  of  2-Hz  Lg,  Sn,  or  coda  to 
predict  mb. 

Variation  is  also  seen  between  stations  consistently  for  all  events  (not  shown).  The 
stations  farthest  onto  the  platform  show  consistently  large  Lg/Pn  and  Sn/Pn  amplitude 
ratios  by  a  factor  of  1. 7-2.5  relative  to  stations  at  high  elevation,  as  well  as  the  largest 
absolute  amplitudes.  The  large  amplitudes  can  be  attributed  to  low  velocities  of  Neogene 
near-surface  sediments  that  lie  beneath  the  platform  stations. 

Conclusions  and  Recommendations 

I.  We  have  shown  that  spectral  fall-off  of  seismic  body  waves  at  2-15  Hz  can 
constrain  attenuation  through  the  Greater  Caucasus  and  through  the  Russian  Platform. 
Overall,  signals  are  2-3  times  more  attenuated  beneath  the  mountains  relative  to  the  adjacent 
shield.  Three-dimensional  inversions  of  t  measurements  confirm  these  observations  and 
demonstrate  that  the  boundary  between  high-Q  and  low-(2  regions  is  sharp  and  closely 
follows  the  northern  slope  of  the  mountain  front.  The  location  of  this  inferred  boundary  is 
somewhat  different  from  that  of  surface  deformation  front,  which  lies  on  the  southern  slope 
of  the  mountains.  Comparison  of  body  wave  with  coda  behavior  suggests  that  intrinsic 


attenuation,  not  scattering,  is  responsible  for  most  of  the  difference  between  mountain  and 
shield  paths.  Scant  laboratory  observations  on  Q  -temperature  relations  imply  70-400"  C 
higher  temperatures  at  the  base  of  the  crust  beneath  the  mountains  relative  to  the  platform. 
These  temperature  increases  are  not  obviously  high  enough  to  trigger  partial  melting, 
although  some  localized  melting  beneath  the  Elbruz  volcanic  region  must  occur.  On  the 
other  hand,  if  the  inference  of  elevated  temperatures  is  correct,  then  a  significant  fraction  of 
the  topography  of  the  Greater  Caucasus  may  be  supported  by  decreased  density  rather  than 
by  crustal  thickening. 

Regardless  of  the  mechanism,  our  results  show  that  attenuation  variations  associated 
with  Alpine-Himalayan  mountain  belts  are  large  and  can  be  abrupt.  Any  attempts  to  use 
amplitudes  of  regional  phases  propagating  near  and  through  these  and  similar  mountains, 
such  as  for  treaty  verification  purposes,  need  to  trace  raypaths  carefully  in  order  to  properly 
account  for  the  sharpness  of  these  boundaries.  The  Q-^  values  derived  here  should  be 
useful  in  calculating  or  correcting  amplitudes  of  seismic  signals  traversing  this  region. 

II.  Extreme  attenuation  beneath  the  Iranian  Plateau  limit  the  utility  of  high-frequency 
waveforms  at  ranges  of  5-15°,  for  sources  smaller  than  M  -  4.5.  Phases  such  as  Lg  are 
absent,  and  the  S  wave  train  is  only  apparent  at  frequencies  of  0.5  Hz.  Reliable  detection 
of  signals  is  likely  to  require  a  combination  of  P  waves  near  1  Hz,  Sn  near  0.25-0.5  Hz, 
and  longer-period  surface  waves.  Discrimants  are  possible  based  on  Sn/Pn  ratios,  but 
likely  Sn  measurements  need  be  made  at  frequencies  of  0.2  -  0.5  Hz  while  Pn  can  be  made 
at  0.5  -  2  Hz.  No  signals  were  detected  at  frequencies  higher  than  4.0  Hz  for  paths  >  10  , 
from  Zagros  events  to  Caucasus  (or  Kopet  Dagh)  stations.  Future  work  should  include  the 
analysis  of  longer  period  surface  waves  as  a  signal  discriminant  along  such  paths,  not 
possible  with  the  instrumentation  in  place  for  the  CNET  experiment. 

III.  Reduction  in  high-frequency  seismic  wave  amplitudes  across  Eurasian  collision 
zones  appears  to  be  associated  with  a  highly  localized  boundary  that  follows  the  northern 
margin  of  the  Greater  Caucasus.  Events  in  the  Racha  region,  on  the  southern  flanks  of  the 
Greater  Caucasus,  show  reduction  in  Lg/P  and  S/P  amplitude  ratios  comparable  to  more 
distant  events  that  cross  the  entire  region.  Although  amplitude  ratios  seem  stable  for  paths 
from  individual  regions,  as  evidenced  by  the  good  correlations  with  mb,  the  attenuation 
effects  are  highly  variable.  A  spatially-varying  attenuation  model  for  the  region  is 
necessary  in  order  to  utilize  Lg/Pn  and  Sn/Pn  discriminants  as  these  amplitude  ratios  vary 
by  over  a  factor  of  5-20  between  different  paths.  The  amplitude  reduction,  seen  throughout 
the  0.5-4  Hz  band,  suggests  that  signals  will  be  small  for  small  regional  events  in  the 
Caucasus-Iranian  Plateau  region,  relative  to  comparable  paths  along  shield  sites.  Robust 
measurements  seem  to  require  both  3-component  analysis  and  averaging  over  several 
nearby  stations,  both  done  here.  Future  work  should  include  quantifying  frequency- 
dependent  0  in  a  way  that  can  be  used  for  path  calibration,  and  adding  complementary 
observations  from  other  networks  and  arrays. 
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Figure  Captions 

Figure  1.  The  Caucasus,  (top)  Regional  topography  in  gray  shades,  500  m  contour 
interval,  network  stations  (solid  triangles)  along  with  telemetry  links,  and  earthquakes  used 
in  Q  measurements.  Plusses  denote  events  for  platform  paths;  solid  circles  denote  events 
for  mountain  paths.  Other  symbols:  stars.  Quaternary  volcanoes  (Philip  et  al.  1989); 
toothed  lines,  major  thrusting  of  the  Greater  Caucasus  (Philip  et  al.  1989;  Triep  et  al. 

1995);  dashed  lines,  location  of  cross  sections  AA’  and  BB’  (Figure  7  and  8).  Inset  shows 
location  of  study  area,  (bottom)  N-S  cross-section  along  AA’,  shown  top.  Also  shown 
location  of  the  network  (CNET),  major  thrusting  and  sediments  beneath  the  stations.  E  = 
volcanic  center  of  Mt.  Elbmz  volcano;  D  =  Dzhirula  massif. 

Figure  2.  Vertical  component  seismograms  recorded  at  station  KTV ,  comparing  regional 
phases  between  mountain  and  platform  paths.  Records  are  plotted  in  same  scale.  Numbers 
in  parenthesis  indicate  group  velocities  in  km/s.  Although  regional  phases  are  efficient  for 
the  platform  path,  they  are  much  weaker  and  longer  period  for  the  mountain  path.  Severe 
attenuation  of  high-frequency  signals  through  the  mountain  paths  may  systematically  bias 
body-wave  magnitude  estimates  here.  Event-distance  (A)  is  approximately  equal. 

Figure  3.  Example  showing  fitting  spectra  procedure  to  measure  t  for  an  event  recorded  at 
station  GUM.  (top)  Seismograms  from  mountain  path  showing  P  and  S  windows  along 
with  pre-event  noise,  barred,  (bottom)  Spectra  for  P  and  S  windows,  labeled.  Inverted 
triangle  shows  estimated /c,  and  /*  values  are  shown  in  box.  Main  arrivals  are  solid  lines, 
and  noise  curves  are  dashed. 

Figure  4.  (left)  r*  estimates  plotted  against  distance  for  mountain  (o),  "mnt",  and  platform 
(+),  "plat",  paths.  Linear  regression  on  r*  (r)  for  events  >0.25°  gives  the  path-averaged  Q 
values  for  each  region,  labeled,  (right)  Variation  in  t  fitting  of  S  for  each  station.  Intercept 
r*  indicates  near-receiver  effect.  Near-site  attenuation  is  systematically  higher  for  MIC, 
NAG  and  KUB,  stations  which  lie  on  soft  sediments.  Station  KNG  has  insufficient 
observations. 

Figure  5.  Variance  (a)  reduction  as  a  function  of  increasing  parameterization  complexity  for 
both  P  (dashed)  and  5  (solid)  waves.  The  parameterizations  are  described  in  Table  2.  Note 
different  scales  for  P  and  S  variances  on  either  side  of  plot.  The  complexity  of 
parameterization,  essentially  the  number  free  parameters,  increases  to  right.  The 
improvement  of  model  parameterization  is  achieved  by  adding  a  low-velocity  layer  beneath 
each  station,  and  by  allowing  Q~'  to  vary  between  mountain  and  platform  regions  (see  text 
for  details). 

Figure  6.  Results  of  attenuation  inversions,  (a)  for  parameterizations  correspond  to 
inversion  H  and  (b)  for  inversion  J.  Each  parameterization  also  shows  number  of  free 
parameters  (F.P.)  that  are  solved  in  the  inversion.  Each  block  is  labeled  by  Q  ‘  values, 
followed  by  formal  errors  (both  multiplied  by  1000).  The  depth  range  of  each  block  is 
indicated  on  the  top  right  of  each  panel.  Triangles  show  station  locations,  and  open  squares 
are  the  Quaternary  volcanic  centers.  Events  used  are  shown  in  Figure  1.  (2  '  foi”  the 
topmost  layer  are  10-100  times  larger  than  those  for  crust  and  upper  mantle  (gray  scales  are 
not  shown).  for  blocks  deeper  than  200  m  is  roughly  equal  to  <2, '  in  each  region. 


Boundary  between  high-Q  and  low-(2  is  sharp  and  closely  follows  the  northern  flank  of  the 
mountain. 

Figure  7.  Cross-sections  (located  on  Figure  1)  for  parameterization  H  and  J.  Format  similar 
to  Figure  6.  Actual  topography  is  shown  at  the  top  of  each  cross-section,  with  lOx  vertical 
exaggeration.  Note  high  <2"'  beneath  mountains. 

Figure  8.  Cross-section  showing  tests  of  resolution  for  Q~‘ .  (top  left)  Initial  model  used  to 
generate  artificial  r*  measurements,  corresponding  to  a  single  lateral  Q  step  that  varies  with 
depth.  Artificial  data  are  generated  for  all  events  used  in  actual  inversions.  Normally 
distributed  random  noise  is  added,  with  variances  equal  to  data  uncertainties  (calculated 
while  determining  t’ ).  Parameterizations  and  a  priori  uncertainties  are  similar  to  actual 
inversions,  (top,  right)  Results  for  parameterization  ff  and  (bottom)  for  J.  Format  similar  to 
Figure  6.  <2~‘  values  are  similar  to  Q~‘  (not  shown).  Original  Q  '  is  recovered  to  within 
18%  for  crustal  blocks  and  54%  for  mantle  blocks,  although  the  overall  correct  pattern 
remains  same.  In  all  cases,  recovered  Q  '  is  less  than  1-(T  formal  error  from  the  correct 
value. 

Figure  9.  Test  for  the  location  of  the  boundary  between  high-Q  and  low-Q  regions,  (a) 
Topography  along  N20°E-S20°W  profile  and  reduction  in  residual  variances  (CJ)  for  t  p  ^ 
t*s  (labeled).  Parameterization  is  similar  to  F  (b),  except  the  location  of  boundary  is  a  free 
parameter.  Inversion  procedure  is  explained  in  text.  For  both  P  and  S,  the  minimum 
residual  crosses  43.265°N  at  43.67“E  (vertical  dashed  lines)  along  the  northern  foothills  of 
the  mountain.  Double  arrows  on  (a)  show  95%  confidence  limits  calculated  from  F-test 
statistics  on  residual  variances.  The  corresponding  attenuation  contrast  between  platform 
and  mountain  is  shown  on  (c). 

Figure  10.  Vertical  records  for  a  regional  distant  mountain  event,  showing  that  the  selected 
time  windows  for  t'  measurements  (horizontal  bars)  are  dominated  by  crustal  signals  {P g 
and  Pt>).  Hence,  spectral  measurements  are  also  likely  to  be  dominated  by  crustal  phases.  A 
Caucasus  velocity  model  (Harvey  et  al.  1994)  is  used  to  identify  Pn,  Pb  Pg  phases 
(dashed  lines). 

Figure  11.  Spectral  measurements  comparing  direct  and  coda  waves,  for  one  seismogram 
traveling  a  shield  path  (top,  A=1.4’,  back-azimuth=165'’)  and  one  traveling  a  mountain  path 
(bottom,  A=1 .5*,  back-azimuth=358*).  The  spectra  of  S  and  5-coda  are  calculated  from 
transverse  component  seismogram,  and  P  spectra  are  from  vertical  component  records. 
Thick  solid  lines  are  the  main  arrivals;  thin  solid  lines  are  the  immediate  coda,  starting  10  s 
after  onset  of  the  main  arrivals.  Dotted  lines  on  P  panels  (left)  are  pre-event  noise;  on  5 
panels  (right)  dashed  lines  are  coda  spectra  starting  at  twice  the  5  arrival  time  (far  coda). 
Triangles  show  the  calculated  source  comer  frequency  for  main  arrivals.  Note  sharper  fall- 
off  of  amplitudes  for  signals  that  traverse  mountains  than  shield.  Corresponding  coda 
spectra  also  have  the  similar  fall-off  characteristics. 

Figure  12.  Map  showing  events  from  PDE  catalog  in  Zagros,  during  period  of  continuous 
CNET  recording.  Open  circles  denote  events  that  did  not  trigger,  closed  circles  show 
events  that  did.  Also  shown  is  path  to  the  one  (of  the  20)  PDE  events  that  triggered  the 
Geyokcha  array. 


Figure  13.  Examples  of  seismograms  from  2  events  in  Zagros,  recorded  at  one  CNET 
station  (NAG).  The  E-W  component  is  shown  to  enhance  S  waves.  Each  panel  shows  the 
raw  trace  (bottom)  and  traces  filtered  in  four  narrow  bands  with  center  frequencies  (from 
bottom  to  top)  of  0.25,  0.5, 1.0,  and  2.0  Hz.  Note  lack  of  Lg  at  all  frequencies,  and 
absence  of  Sn  at  frequencies  below  0.5  Hz. 

Figure  14.  Magnitude  versus  source-receiver  range  for  the  20  Zagros  events  in  PDE 
catalog  during  the  time  QIET  continuous  records  were  collected.  Each  event  is  denoted  by 
2-3  circles,  showing  closest  and  farthest  network  station.  Solid  symbols  showed 
detectable  signal  or  resulted  in  an  event  trigger,  while  open  symbols  did  not. 

Figure  15.  Records  from  two  events  from  ~3°  distance  recorded  at  CNET.  Waves  for  the 
upper  event  in  Dagestan  travel  largely  on  the  Russian  Platform  while  waves  for  the  Lesser 
Caucasus  event,  bottom,  traverse  the  Greater  Caucasus  and  Kura  Basin.  Seismograms  are 
labeled  by  arrivals  of  Pn,  Sn,  and  group  velocities  in  km/s.  Notice  absence  of  late-phase 
energy  on  lower  traces,  particularly  in  and  after  Lg  window  (3.6  km/s).  The  upper  event 
shows  small  Pn  relative  to  Sn  and  late  phases,  and  slow  high-frequency  surface  waves  (<3 
km/s).  Map,  lower  right,  shows  network  configuration  and  event  locations. 

Figure  16.  Ratios  of  RMS-Lg  amplitude  to  P  amplitude  for  events  1-10°  from  CNET, 

1992.  Each  panel  represents  a  different  frequency  bandpass,  labeled  top.  Amplitude 
calculation  is  described  in  text.  Symbols  are  centered  on  events  and  represent  network- 
averaged  amplitudes.  Size  of  symbols  scales  to  amplitude  ratios  according  to  key  ~  left- 
hand  and  right-hand  numbers  are  for  left  and  right  column  of  figure,  respectively. 

Figure  17.  Variation  in  RMS  amplitude  ratios  as  a  function  of  frequency  for  three  regions. 
Phases  used  are  ratio  of  Sn/Pn  (top)  and  Lg/Pn  (bottom)  for  each  of  3  subregions.  Path 
lengths  are  3  -  10°.  Regions  on  Figure  18. 

Figure  18.  (left)  Three-component  RMS  amplitudes  of  2-Hz  wave  train  vs  mb.  Coda  is 
15-30  s  window  centered  at  twice  the  S  travel  time.  Lg  window  spans  3.6  -  2.6  km/s  group 
velocity,  and  Pn  window  extends  to  25%  of  S-P  time.  Signal-to-noise  ratios  exceed  2.0 
for  all  measurements,  and  measurements  are  averaged  over  3  or  more  network  stations. 
Different  symbols  correspond  to  regions:  DA  -  Dagestan;  RA  -  Racha  aftershock  zone;  TU 
-  Turkey.  The  m/,  value  is  taken  from  PDE  for  TU  and  JSPC  measurements  elsewhere 
(Harvey  et  al.,  1994).  Correlation  coefficients  between  log(ainplitude)  and  are  shown 
in  upper  left  of  each  plot.  Amplitudes  are  corrected  for  cylindrical  geometric  spreading  to 
nm/s  at  100  km.  (top  right)  Map  of  regions;  GC  =  Greater  Caucasus,  LC  =  Lesser 
Caucasus,  FSU  =  Former  Soviet  Union. 
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Table  2  :  Description  of  inversion  parameterization 

Variances.  X 10  s 


Description _ , 

homogeneous  half-space 
2  layer,  1  block  in  each  laya 

1  layer  and  2  blocks 

2  layer,  2  blocks  in  each  layer 
2  layer,  1  block/station  in  top 
homogeneous  lower  layer 

2  layer,  1  block/station  in  top 

2  bloc^  in  lower  layer 

3  layer,  1  block/station  in  top 

2  cmstal  blocks 

1  mantle  block 

3  layer,  1  block/station  in  top 
lower  layers  similar  to  F 

2  layer,  1  block/station  in  top 
lower  layer  thin  blocks 

3  layer,  1  block/station  in  top 
lower  layers  similar  to  I 

2  layer,  1  block/station  in  top 
80  km  X  80  km  block  lower 


#  free  parameters 

P-wave 

S-wave 

1 

0.705 

1.543 

2 

0.541 

1.474 

2 

0.578 

1.115 

4 

0.463 

1.091 

7 

0.464 

0.971 

8 

0.380 

0.637 

9 

0.381 

0.675 

10 

0.375 

0.635 

17 

0.369 

0.638 

28 

0.341 

0.609 

48 

0.341 

0.618 

Tstole  3  :  Events  from  S.  Iran  found  In  Caucasus  and  Geyotoha  triggered  caaiog 
latitude  longitude  origln.trne  depth  delta 

30.0070  52.24^  1 0/13/1 9M  33.0000  15.880 

(286) 

0:11:04.700 _ 


29  7730  51.1030  10/23/1993  33.0000  1  5.683 

(296) 

13:26:42.100 _ _ 


27  8880  57.4150  11/01/1993  33.0000 

(305) 

_ 6:42:00.800 _ 

28  7590  51.9900  11/03/1993  33.0000  16.920 

(307) 

_ 16:42:17.600 _ 

27  9200  57.2340  3/03/1994  33.W00  19.835 

(062) 

_ 14:56:46.900 _ 

28  9000  52.4650  3A)3/1994  33.0000  16.962 

(062) 

23:54:01.000 _ 


29.0260  52.7260  3A)7/1994  10.X00  16.945 

(066) 

0:57:19.600 _ 


33.1460  48.0330  3A)7/1994  26.0000  11.587 

(066) 

10:54:57.400 _ 


28  6330  53.1290  3/19/1994  33.0000  17.452 

(078) 

4:54:00.300 


28  7270  52.5880 


signals  seen?  compared  to  sta  list  in  triggered 

_  triggered  NET  catalog 


Caucasus 


KIV  <&  1  hz  I  Caucasus  j  NAQ 


QUM  KIV  MIC  Caucasus 
NAQ  KNQ  #  2 
H2 


Caucasus 


28  9940  I  52.7450 


28  87i0  52.5960  3/23/1994  33.0000  17.034 

(082) 

17:14:45.300 _ 


29  096'o  51.2560  3/29/1994  33.0000  16.359 

(088) 

7:56:53.900 _ 


3/30/1994  54.0000 

(089) 

19:55:46  000  _ 


20  9940  52  7450  3/30/1994  54.0000  7  38 

(069) 

19:55  46.000 _ 


28  8210  52  7450  4A)3'1994  23  0000  17  137 

(093) 

6:51:57  500 _ 


20  8870  52  7050  4A)3/1994  33  0000  17  062 

(093) 

7:19:36  300 _ 


27  4770  54  3850  4/21/1994  33  0000  18  979 

(111) 

11:50  32  600 _ 


30  9990  50  6350  4/26/1994  33  0000  14  387 

(116) 

6:26:01,000 _ 


4/30/1994  330000 

(120) 

23:11:13  400 


5A)3/1994  33  0000 

(123) 

20:20:18800 


29  0080  i  52  6870 


33  5010  I  473840 


yes  @  2  Hz 


yes  @  2  Hz 


yes  @  2  Hz 


yes  @  2  Hz 


yes  2  Hz 
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